Identification of H2CCC as a diffuse interstellar band carrier 
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ABSTRACT 

We present strong evidence that the broad, diffuse interstellar bands (DIBs) at 4881 and 
5450 A are caused by the B 1 Bi <- X 1 A i transition of B 2 CCC (7-C 3 H 2 ). The large widths of the 
bands are due to the short lifetime of the B 1 Bi electronic state. The bands are predicted from 
absorption measurements in a neon matrix and observed by cavity ring-down in the gas phase 
and show exact matches to the profiles and wavelengths of the two broad DIBs. The strength 
of the 5450 A DIB leads to a Z-C 3 H 2 column density of ~5xl0 14 cm" 2 towards HD 183143 and 
~2xl0 14 cm~ 2 to HD 206267. Despite similar values of E(B - V), the 4881 and 5450 A DIBs 
in HD 204827 are less than one third their strength in HD 183143, while the column density of 
interstellar C3 is unusually high for HD 204827 but undetectable for HD 183143. This can be un- 
derstood if C3 has been depleted by hydrogenation to species such as Z-C3H2 towards HD 183143. 
There are also three rotationally resolved sets of triplets of Z-C3H2 in the 6150— 6330 A region. 
Simulations, based on the derived spectroscopic constants and convolved with the expected in- 
strumental and interstellar line broadening, show credible coincidences with sharp, weak DIBs 
for the two observable sets of triplets. The region of the third set is too obscured by the a-band 
of telluric 2 . 



Subject headings: ISM: general 
techniques: spectroscopic 

Introduction 



Diffuse absorption features were first noted in 
the spectra of 'red' B stars some 90 years ago, 
but only recognized as arising from the interstel- 
lar medium in the 1930s. Initially, broad bands 
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operated by the National Research Council of Canada, the 
Institut National des Sciences de l'Univers of the Centre 
National de la Recherche Scientifique of France and the 
University of Hawaii, and at the Dominion Astrophysical 
Observatory, Herzberg Institute of Astrophysics, National 
Research Council of Canada. 
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some tens of A wide were identified, but with im- 
proved instrumentation many hundreds were sub- 
sequently found — mostly in the yellow/red with 
each being only a few A wide. Although the 
strength of the diffuse interstellar bands (DIBs) is 
correlated roughly with colour excess, none shows 
the characteristics expected for absorption within 
the dust grains themselves. Rather, they are pre- 
sumed to arise from gaseous molecules. Their 



characteristics have been summarised ( Snow & 



McCall 2006). Several hypotheses have been put 
forward for their origin; most prominent being 
carbon chains and polycyclic hydrocarbon cations 
( Herbig|2000 ). Advances in laboratory techniques 



over the last decade have enabled the measure- 
ment of a number of such electronic spectra in 
the gas phase. These include some bare carbon 



chains (Maier et al. 2004), those containing one 



hydrogen, nitrogen and their ions ( Motylewski 



et al. 2000), and a few aromatic cations (Bicn- 



nier et al. 2003 ) . Individual coincidences of three 



very weak DIBs with the gas-phase absorptions 



of naphthalene cation (Iglcsias-Groth et al. 2008) 



and, recently, for one band of diacetylene cation 
( Krelowski et al.||2010 ) have been reported. 

A number of di- and triatomics are known to 
be present in the diffuse clouds: CH, CH + , CN 



as well as HJ and C 3 (see Snow & McCall 2006 



for references therein). In addition, mm- wave ob- 
servations indicate the presence of more complex 



carbon containing molecules, e.g. 0C3H2 (Cox et 
aLl[T988l ), C 4 H ( [Liszt fc LucaspOOO ). 

Recently, a match between a broad DIB, at 
5450 A with -10 A FWHM, and the absorption 
of an unknown species produced in a plasma dis- 



charge of acetylene was reported (Linnartz et al 



2010). The available experimental data pointed 
to a molecule containing only carbon and hydro- 
gen. Here, we present evidence that the molecule 
responsible is propadienylidene, H2CCC — linear 
three carbon atoms with two hydrogens off-axis 
on the terminal carbon; henceforth and usually 
designated as Z-C3H2. The molecule has two 
other broad absorption bands at around 4890 and 
5170 A, close to known DIBs, as well as rotation- 
ally resolved bands at wavelengths above 6000 A. 
Narrow DIBs of the expected strength are also 
found at the positions of two of the triplet systems 
predicted to be centred near 6159 and 6252 A. The 
spectral region for the third system, weaker than 
the other two, is too confused by strong telluric 
oxygen features. All of this is evidence for an iden- 
tification for a carrier molecule which gives rise to 
a series of DIBs. 

The recognition that the broad band seen in 
the lab which coincided with the 5450 A DIB 
was due to absorption by Z-C3H2 became evident 
from a number of previous studies, followed by 
the ones reported here. In its electronic ground 
state X 1 Ai, Z-C3H2 (C2 V symmetry) has been 
well characterized by mm-wave spectroscopy in 
the laboratory ( Vrtilek et al.||1990 ) and identified 



by this signature in dense interstellar clouds ( Cer- 



(Cernicharo et al. 


1999 


). It is an isomer 


ubiquitous cyclic C3H2 ( 


Madden et al.|1989 



electronic absorption spectrum of Z-C 3 H 2 was first 



observed in a neon matrix at 6K (Hodges et al 



2000 ) . In the visible region it shows two electronic 
systems; a weak dipole forbidden A 1 A2-s—A 1 Ai 
transition which becomes partly allowed by vi- 
bronic effects with the near lying B 1 B 1 electronic 
state. Then it shows an intense _B 1 B 1 <— X 1 A X 
transition in the 4000— 5500 A range. The spec- 
trum has now been re-recorded following mass- 
selected deposition of C3HJ cations followed by 
neutralization; the earlier measured spectrum was 
obtained by photolysis of a chemical precursor in 
the matrix. The spectra are by and large identical, 
except that the present recording (Fig. 1) shows 
better resolution of the bands and there are wave- 
length errors for some of the reported absorptions. 
Due to matrix-gas shifts, a direct comparison with 
DIB data is inconclusive. 

Thus a goal a number of years ago was to ob- 
serve the B 1 Bi 4- X 1 Ai system of Z-C3H2 in the 
gas phase. Three techniques were employed: ab- 
sorption measurements in a supersonic jet expan- 
sion through which a discharge runs and detection 
with cavity ring-down (CRD) either with pulsed- 
or continuous lasers, as well as resonance enhanced 
multi-photon ionisation. The experiments were 



unsuccessful (Achkasova et al. 2006), from which 
it was concluded that the B i B 1 excited electronic 
state has a lifetime less than a picosecond, causing 
individual line broadening of over 10 cm -1 . How- 
ever, at the same time, three rotationally resolved 
bands lying in the 6000-6600 A region were ob- 



served (Achkasova et al. 2006 Birza et al. 2005) 



In the latter two articles, the rotational analyses 
proved that these belong to Z-C3H2, correspond- 
ing to transitions from the v=0 level of the ground 
X : Ai state to vibronic levels in the A X A2 state, 
lying near the origin of the B 1 B 1 electronic state. 
The bands in the dipole forbidden A 1 A 2 4- X 1 Ai 
transition gain intensity by vibronic mixing with 
the B 1 Bi levels. The rotational analysis is un- 
ambiguous. These bands, lying around 6159, 6252 
and 6319 A in the gas phase show a complex struc- 
ture characteristic of an asymmetric top. 



nicharo et al. 1991), and in the diffuse medium 
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2. Laboratory Observations 
2.1. Techniques 

The two laboratory techniques used were ab- 
sorption measurements in a neon matrix and CRD 
in the gas phase. For the former C3HJ was 
produced from allene in an electron impact ion 
source. The mass-selected m/e=38 ions were then 
co-deposited with excess of neon onto a sapphire 



neutralization. That this is the absorption of lin- 
ear C3H2 has been proven (Hodges et al 



2000) 



surface at 6K to form the matrix (Freivogel et 
al. 1994). This took a few hours using nA ion 



currents to produce a matrix of 150 fim thickness, 
spread over a few cm 2 . The ions were subsequently 
neutralized by UV irradiation of the matrix, re- 
leasing electrons from the counter ions present. 
The absorption spectra were measured by passing 
light through the thin side of the matrix over a 
path length of 2 cm, dispersion of the transmitted 
light by a grating and detection with a CCD cam- 
era. Another absorption spectrum was obtained 
using allene as precursor in the matrix followed by 
xenon lamp VUV (1470 A) photolysis. The mass- 
selected measurements on /-C3D2 were achieved 
with dideutero-acetylene in the ion source. 

In the gas phase Z-C3H2 was produced in a su- 
personic free-jet expansion discharge. As precur- 
sor, 0.5—1% allene in argon, or acetylene and di- 
acetylene in helium, was used. The CRD set up 
is as described (Linnartz et al. 1998). The dis- 



charge was operated at 10 Hz while the laser was 
fired every 20 Hz. The CRD signal, which is not 
background-free, could then be corrected for base- 
line. It was found that the strongest signals were 
observed with allene, but at the same time the 
vibrational and rotational temperatures attained 
were a bit higher than with diacetylene. Further- 
more, allene produced other species causing broad 
absorptions in the range of interest (e.g. I-C3H). 
Dideutero-acetylene in helium was used in the dis- 
charge expansion for the measurements on Z-C3D2. 

2.2. Laboratory Results 

We focused on whether the 5450 A absorption 
band coinciding with the DIB at this wavelength 
( |Linnartz et al.|2010| ) could be the B 1 B i <- X 1 Ai 
transition of I-C3R2, he. corresponding to the 
most intense band in the 6 K neon matrix spec- 
trum at 5417 A (Fig. 1). The middle trace was 
obtained using C3HJ mass selection followed by 



Our spectrum is the same, but with a better reso- 
lution; for example the site-structured bands show 
two components and the correct calibration of the 
5445/5417 A band is given in Table 1 . The spec- 
trum below 4700 A is not shown, because overlap- 



ping bands are present there due to Z-C3H (Ding 
et al.||200~ ), produced as a result of H loss during 
the deposition of C3HJ. 

The top trace is the spectrum obtained via a 
traditional matrix approach: allene was embed- 
ded in neon and following photolysis using a low- 
pressure xenon lamp the system appears. This 



spectrum is identical to the one published ( Hodges 
et al.|2000 ), though produced by a different chem- 
ical scheme and irradiation. It is included because 
the intensity ratio of the respective bands is more 
reliable than in the mass-selected spectrum where 
the absorption involves much smaller concentra- 
tions and thus background correction is not so re- 
liable. The doublet at 5445/5417 A is not resolved 
but is indicated from the asymmetry of the band. 
The structure of the 5417 A band is due to matrix 
site effect, because the same asymmetric pattern is 
seen for the other bands (5143, 4856 and 4633 A). 

Bottom trace in Fig. 1 shows the spectrum mea- 
sured after deposition and neutralization of C^D^ 
ions into a neon matrix. The site structure is not 
present for the prominent band at 5412 A; such an 
effect is often seen on isotopic substitution in ma- 
trices. The band at 5412 A is shifted by 17 crn^ 1 
towards higher energy in comparison to the 5417 A 
band of Z-C3H2. 

The assignment of the bands to vibrational ex- 
citation in the upper electronic state of Z-C3H2 is 
established (Hodges et al. 2000). All the transi- 



tions originate from the v=Q level of the X : Ai 
ground state because the temperature of the neon 
matrix is 6K. The strongest 5417 A band corre- 
sponds to the excitation of the totally symmetric 
V2 (C=C stretching) mode in the B 1 Bi state (2j 
transition). Other strong bands in the spectrum 
of C3H2 are due to its overtone (2g at 4856 A) and 
combination with the V4, mode (2j4j at 5143 A) 
as indicated in Fig. 1. The origin band of the 
B 1 Bi<— X X K\ transition is weaker and is influ- 
enced by vibronic interaction of B l YS\ with the 
nearby lying A 1 K-i state. It is located somewhere 
within the group of weak bands above 6000 A (be- 
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tween a and c in Fig. 1). This has been considered 



theoretically (Hodges et al. 2000) and also dis- 



cussed in the analysis of the rotationally resolved 



spectra in the gas phase in this region ( Achkasova 
eFaLl[2006l ). 



CRD is a sensitive detection method; however, 
the lack of mass selection is a drawback. There- 
fore, the matrix spectra of the mass-selected C3H2 
and C3D2 serve for the identification of the bands 
in the gas-phase CRD experiments. Though the 
wavelength shift going from neon matrix to the 
gas phase can be 10— 50 A for electronic transi- 
tions, the separation of vibrational peaks ra rely 



differ by more than a few A within one state (Ja- 



cox 



1994). Thus, by knowing the wavelength of 



the 2q transition in the gas phase, the positions 
of other bands (e.g. 2q or 2j 4q) can be predicted 
from the separation of the peaks in the neon spec- 
trum (Fig. 1). However, this procedure cannot be 
applied for the 2q band of Z-C3H2 itself because 
the separation from the origin to the 2q peak is 
not known. 

The aim was therefore to detect the gas-phase 
band of Z-C3H2 that corresponds to the 2j transi- 
tion in the neon matrix at 5417 A. The known ro- 
tationally resolved bands of Z-C3H2 above 6000 A 
were first detected using acetylene, diacetylene 
or allene seeded in argon in the discharge. The 
strongest signal was observed with allene. The 
spectrum was then measured in the 5420±50 A 
range as indicated by the 2g absorption in neon. 
Broad absorptions were detected with all the pre- 
cursors and under the same discharge conditions 
as for the rotationally resolved A > 6000 A bands 
of Z-C3H2. The spectrum obtained with allene pre- 
cursor is shown in Fig. 2 a. It consists of two broad 
absorptions, whose relative intensities vary with 
the precursor proving that they originate from two 
different molecules. The intensity of the band at 
lower wavelength (—5450 A) correlates with the 
identified absorptions of I-C3H2 in the range above 
6000 A, but the one around 5480 A does not. 

The absorptions in Fig. 2 a have been decon- 
volved using the two Gaussians shown. The 
central wavelength of the first peak is —5450 A 
(FWHM —13— 17 A), whereas the other is centered 
at 5473-5476 A (FWHM -32-42 A). The promi- 
nent C2 absorption lines have been partially re- 
moved from the spectrum by taking a wavelength 
window of ten points, finding the minimum and 



then repeating the procedure at the following po- 
sition of the spectrum. By this means the broad 
underlying absorption is retained, and fitted by a 
Gaussian. Fig. 2 a shows the superposition of the 
two Gaussians and the quality of the fit. 

Gas-phase measurements on Z-C3D2 were also 
carried out. In the spectral window of inter- 
est only one broad band (central wavelength at 
5458(3) A and FWHM -19-24 A) has been de- 
tected using dideutero-acetylene as the precursor 
(Fig. 2 b). The FWHM does not differ much from 
that of the 5450 A band. The position of this band 
is shifted —27 cm -1 to the lower energy relative 
to the 5450 A band, similar to the neon matrix 
change (17 cm -1 ). 

Comparison of the neon matrix spectra of l- 
C3H2 and Z-C 3 D 2 (Fig. 1) with the CRD spectra 
(Fig. 2) leads to the conclusion that the 5417 A 
band in the neon matrix and that at 5450 A in 
the gas phase originate from the same molecule, 
namely Z-C3H2. Also the fact that the intensity 
of the band at 5450 A observed in the gas phase, 
correlates with the intensity of other, already iden- 
tified bands of Z-C3H2 in the rang e above 6000 A, 
substantiates this conclusion. 

From the position of the 2j band of Z-C3H2 in 
the gas phase (5450 A) the wavelengths of the 2q 
and 2q 4j transitions in the gas phase can be pre- 
dicted using the separation of the peaks in the 
neon spectrum. For this the most intense peaks 
are taken (5417, 5143 and 4856 A), because the 
lowest energy site is not resolved (though its pres- 
ence is indicated by the asymmetry of the peak) . 
The difference between the 2j and 2q band max- 
ima is 2130 cm -1 and this leads to a gas phase 
prediction of 4883 A for the 2q transition; the un- 
certainty is difficult to guess because of broaden- 
ing mechanisms in the matrix. The 2j 4j band is 
located 984 cm -1 above the transition in the neon 
matrix and in the gas phase it is then in the 5173 A 
wavelength region. 

In the region around 4883 A predicted for the 
2q band, broad absorptions are detected (Fig. 2 c) 
with narrower components imposed on them. The 
four bands at 4913, 4870, 4856 and 4844 A belong 
to Z-C3H because they have been identified in ear- 



lier mass-selected studies in the gas phase (Ding 
et al. 2001). After subtraction of the Z-C3H ab- 



sorptions (the 4 Gaussians shown in Fig. 2 c) from 
the laboratory recording, the residual broad band 
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is obtained (Fig. 2 d) . The Gaussian fit is centered 
at 4887(3) A with FWHM of -25 A. Its position 
matches the 4883 A wavelength predicted from the 
neon matrix spectrum of Z-C3H2. The FWHM of 
this band is somewhat larger than that of the 2j 
band at 5450 A probably due to a shorter lifetime 
in the V2—2 vibrational state of B 1 Bi. These ar- 
guments, together with the fact that the 5450 and 
4887 A bands were observed under the same ex- 
perimental conditions as the rotationally resolved 
transitions of Z-C3H2 above 6000 A provide the ev- 
idence that these absorptions originate from the 
Z-C3H2 molecule. 

The detection of the band around 5173 A, cor- 
responding to the 2j 4q transition of Z-C3H2, was 
not successful. This is because this transition is a 
factor of three weaker than 2q and its intensity is 
below the detection limit of the CRD set-up. 



It was already established (Linnartz et al. 2010 ) 
that the 5450 A gas-phase absorption coincides 
with a DIB at this wavelength. This band is now 
assigned to the 2j B «— X 1 Ai transition of l- 
C3H2. A DIB corresponding to the 2g band in the 
gas-phase observations is present (see Section 4). 

3. Simulation of Band Profiles 
3.1. Bands Above 6000 A 

Three systems in the 6150— 6330 A region have 
been subject to two earlier high-resolution, gas- 



phase studies (Achkasova et al. 2006 Birza et al 



2005 ) . The rotational analysis is unambiguous and 



proves that the absorptions belong to Z-C3H2. 

Z-C3H2 has a near prolate top structure with 
C%v symmetry. The rotational constants in the 
A x Ai ground state are ,4=9.6451, 5=0.3533 
and C=0.3404cm- 1 . The K=l - K =0 and 
K=2 — K=0 separations are ~10 and ~40 cm -1 
respectively. Thus at very low temperatures only 
the K=0,1 levels are populated. According to the 
selection rules only transitions with AAT=±1 are 
observed. 

In Fig. 3, left, are shown the simulations of the 
bands of Z-C 3 H 2 above 6000 A, using the derived 
spectroscopic constants and a Lorentzian profile 
of 0.15 A width, corresponding to the resolving 
power 40,000 used in astronomical measurements, 
for each individual line. The pattern is shown 
for temperatures of 5, 20 and 60 K. Three com- 



ponents, K=2^- 1, K=l^Q and K=0^l would 
be observed. At 5 K the population becomes con- 
centrated in the K=0 level of the A 1 Ai electronic 
ground state and thus the middle peak dominates. 
As the temperature rises the outer two peaks gain 
intensity because the population of the K—\ level 
increases. 

Though the resolving power in many astro- 
nomical measurements is around 40,000, observa- 
tions such as toward HD 206267 sample several 
clouds as the profiles of K and Ca lines show 



(Pan et al. 2004[ ) and this leads to broadenings 
of 10-20 km s _i , i.e. 0.2-0.4 A in this wavelength 
region. Thus the absorption pattern of Z-C3H2 in 
such clouds would look like that shown in Fig. 3, 
right, where a 0.4 A width (Lorentzian) is used in 
the simulation. 

3.2. 5450 A Band 

The 5450 and 4887 A absorption bands of l- 
C3H2 are broad, without any rotational struc- 
ture. This arises due to the short lifetime of 
the B x Bi electronic state as a result of fast in- 
tramolecular relaxation process to the nearby ly- 
ing A 1 A% state. The profiles of the broad 5450 
and 4887 A bands can be predicted by simulating 
the 2j and 2q B 1 Bi <- X 1 Ai transitions with the 
known rotational constants fo r the bands above 
6000 A ( |Achkasova et al.|2006[ ). Due to the broad- 
ness of the features, the precise value of the rota- 
tional constants has little influence on the overall 
contour. The bottom trace in Fig. 4 shows the case 
with a 1 A linewidth corresponding to 1 ps lifetime 
and temperature of 10 K. Three components in the 
profile of this band would still be seen. 

The best fit to the observed 5450 A profile was 
obtained with a 10 A linewidth and 60 K tempera- 
ture (top trace in Fig. 4) . The temperature of 60 K 
is chosen because the rotationally resolved bands 
of Z-C3H2 measured above 6000 A fit to 50-70 K. 
The linewidth of 10 A corresponds to a 100 fs life- 
time, which is often encountered in the excited 
electronic states of polyatomic molecules subject 
to conical intersections. In fact in the case of l- 
C3H2 the next higher lying C : Ai electronic state, 
with onset near 2500 A, has been shown experi- 



mentally to have a lifetime of 70 fs (Noller et al 



2009), so the inferred value of 100 fs is not unrea- 



sonable. 
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4. Astronomical Observations 



4.2. Narrow Triplets Above 6000 A 



4.1. Broad Bands 

Broad absorptions are observed in the lab at 
4887 and 5450 A in the gas phase and one pre- 
dicted in the 5165-5185 A range but expected to 
be only one half to one third the strength of the 
5450 A absorption. The absorption spectrum in 
Fig. 1 is a linear measurement and thus predicts 
that any corresponding DIBs should have compa- 
rable EWs, and FWHM. There is a strong, broad 
DIB at 4881 A with twice the FWHM (~25A), 
but similar central depth to the 5450 A DIB. It 
is listed in many DIB compilati ons (| Jenniskens "fc| 
Desert|[l994l |Hobbs et al.|[2008| [Herbig||1995| ) and 
was recorded anew as part of this study. Unfortu- 
nately, stellar H/3 (4861 A) contaminates its short- 
wavelength wing. 

Spectra were acquired in the regions of these 
DIBs at the Dominion Astrophysical Observatory 
with the 1.2 and 1.8 m telescopes at a resolution 
of i?=18,000 in early 2010. The reddened target 
stars are listed in Table 2 and the spectra displayed 
in Figs. 5 and 6 for the 4881 and 5450 A regions, 
respectively, together with those of unreddened 
stars. The spectra for individual stars were aligned 
to the mean rest wavelength of the interstellar K I 
line in each case. 

The coincidence of the lab and DIB profiles in 
the cas e of the 5450 A band has already been dis- 
cussed ( Linnartz et al.|2010 ). The Gaussian (25 A 
FWHM) obtained from the CRD lab measurement 
of the 4887 A absorption (Section 2.2) is superim- 
posed on the DIB spectra in Fig.[5j The agreement 
is again excellent. The wings of the strong stel- 
lar H/3 line at 4861 A makes that part of the DIB 
hard to define for HD 46711. For the supergiant, 
HD 183143, there might also be a broad emission 
component. 

The case for a FWHM~20A DIB in the 
5165—5185 A range is complicated by the presence 
of a previously known, and well defined, broad 
DIB centred at 5160 A. Unfortunately, the pres- 
ence of strong stellar lines near 5173 A in all of 
our target stars made it impossible to say any- 
thing about an interstellar absorption feature in 
that region. 



The triplets offer quite a different challenge 
from the broad bands when it comes to unambigu- 
ous detection. An exact wavelength match was 
less critical for the broad bands. There are many 
fewer broad DIBs known than sharp ones mak- 
ing chance coincidence less likely. Also, the broad 
DIBs can be matched in width without significant 
distortion from spectral resolution and blurring by 
the velocity dispersion in the interstellar clouds. 
For the 4881 and 5450 A regions there is little 
contamination by telluric lines and stellar lines 
are easily distinguished. By contrast, the spec- 
tra of the late-type B giants are full of stellar lines 
which are hard to distinguish from sharp DIBs. 
The highest density of sharp DIBs is in the 6000 A 
region, one every few A, which makes the chance 
of coincidence likely. There are also many telluric 
lines. It should also be noted that the sharp DIBs 
of interest have central intensities of <1 %. 

As a result of these obstacles, we were restricted 
to reddened early B or O stars and spectra with 
a resolution >40,000 and S/N of several thousand 
per A. We drew on CFHT archival spectra and 
others were especially acquired of HD 179406 with 
the CFHT ESPaDONs spectrograph in June 2010. 
The stars are listed in Table 3 where 'Gee' cor- 
responds to Gecko with a median resolution of 
120,000 and 'Esp' is ESPaDONs with R -80,000. 
Both are bench spectrographs at the CFHT. It 
should be emphasized that the archival spectra 
were optimized for other programs. Low order 
polynomials have been applied to remove residual 
undulations in the final spectra. 

All of the target stars are seen through mul- 
tiple interstellar clouds with some differential ve- 
locity ranges >10kms~ 1 as can be seen in Fig. [7] 
where the interstellar KI line profiles are plotted 
and aligned on their strongest components. This 
alignment was applied to the individual stellar 
spectra and the corresponding interstellar veloc- 
ities are listed in Table 3. 

Figs. [8] and [9] display the spectra of the reddened 
stars for two sets of the predicted Z-C3H2 triplets 
together with model absorption profiles for 10 K 
and a 0.15 A FWHM, the latter to account for in- 
strumental and interstellar broadening. There are 
clearly DIBs within the limits of the model profiles 
at 6159, 6166, 6244, 6251 and 6259 A. The line pre- 
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dieted at 6152 A was expected to be significantly 
weaker than those at 6159 and 6166 A so the ab- 
sence of a DIB at that wavelength is not necessar- 
ily a contradiction. The DIBs at 6244 and 6251 A 
appear to be double with the stronger component 
corresponding quite closely to the predicted ab- 
sorption maximum. The comparison between the 
lab and interstellar wavelengths together with the 
EW of the individual astronomical DIB for each 
of the target stars are given in Table 4. 

For the triplets predicted near 6320 A, the re- 
gion is heavily obscured by the telluric a-band of 
O2 which, coupled with velocity corrections to the 
barycentre and alignment to the interstellar KI 
velocities, makes it impossible to detect any DIB 
which might coincide with the predictions. 

In none of the coincidences do the predicted 
profiles correspond well to those of the DIBs, sug- 
gesting that the intrinsic broadening is probably 
greater than the 0.15 A assumed for the model. 
While the coincidences between DIBs and the 
model profiles are not unequivocal given the high 
density of sharp DIBs in this spectral region, their 
absence would have cast some doubt on the case 
Z-C3H2 as a DIB carrier. 

5. Discussion and Conclusions 

The choice of a rotational temperature for l- 
C 3 H 2 in diffuse clouds has a consequence for the 
comparison with the laboratory data, as the sim- 
ulations in Figs. 3 and 4 show. I-C3H2 has a large 
dipole moment (around 4D) and may cool effi- 
ciently, in much the same way as do the diatomics 
CH, CH+ and CN, which attain 2.7 K in the dif- 
fuse medium. On the other hand, the tempera- 
ture inferred for other polar polyatomics detected 
in the diffuse medium are often higher, for exam- 
ple arou nd 50 K for amm onia by mm-wave spec- 



troscopy ( Liszt et al. 2006 ) and H3 by IR measure- 



ments ( McCall et al. 2002 ) . Due to the broadness 
of the 5450 A DIB the interstellar temperature of 
Z-C3H2 can not be clearly determined; anything in 
the range 10— 60 K would fit as the two simulated 
profiles with 10 A linewidth in Fig. 4 show. 

As a follow up to the identification of ^-C3H 2 
as a DIB carrier, one can make an estimate for 
the column density on the basis of the broad 
5450 A DIB. There are two calculations (Mebel 



of the i? 1 Bi<— X x Ai electronic transition, 0.016 
and 0.009. The former value was obtained with a 
higher level method and is used for the evaluation. 

The / value for the 2j transition (5450 A band) 
within this system is estimated from the calcu- 
lation of the Franck— Condon factors; these are 
shown in Fig. 5 of ( |Mebel et aL|l998| ). This gives 
/(2o)//totai(=0.016) = 0.2. Using the expression: 



iV(cm 



= 1.13 x 10 



, EW(A) 
A(A) 2 x / 



yields for A=5450A, /(2j)=0.003, EW=0.4A to- 
wards HD 183143: 7V(Z-C 3 H 2 )«5x 10 14 cur 2 . 
Towards HD 183143 the column den sity of Hj|~ 



has been determined to be 10 14 cm 2 (McCall et 
al.|2002| ) and 7V(C 2 )«10 14 cm~ 2 ( |Oka et al |2003f r 
Thus the N value estimated for Z-C3H2 is not un- 
reasonable. 

An obvious question is the role of c-C3H 2 . Ra- 
dioastronomy has established that this cyclic iso- 
mer is ubiquitous in dense and diffuse interstel- 
lar media. But C-C3H2 does not have an allowed 
electronic transition in the DIB region. Its ab- 
sorption has been observed in the labora tory, in 
a neon matrix with onset around 2700 A (Seburg 



et al. 1997). Apart from Z-C 3 H 2 , also Z-C3H has 



many absorption bands in the DIB region (Ding 



et al. 2001 1. Some of the peaks are actually ob- 



served in the CRD spectra measured, nearby to 
the broad bands of Z-C3H2. However, as has al- 
ready been pointed out, due to the small oscil- 
lator strength of such individual transitions the 
DIB absorptions would have EW of mA or less. 
Nevertheless, chemically the presence of I-C3H in 
addition to the here proven I-C3H2 could be inves- 
tigated with high S /N DIB measurements because 
the gas-phase data are available. Another related 
aspect concerns an earlier discussion on the pos- 
sible presence of I-C3H2 in the diffuse medium 



(Guethe et al. 2001). It was pointed out that l- 



et al. 1998) of the overall oscillator strength, /, 



CsH 2 formation could be significant via associate 
electron attachment to Z-C3H2. 

As the 4881 and 5450 A DIBs correspond to the 
two broad absorptions of ^-C3H2, a further aspect 
has to be considered, namely the estimate of the 
expected EW for the narrow b, c bands around 
6251 A and 6159 A with the triplet structure in 
the DIB observations as discussed in Section 4.2. 
The matrix absorption spectrum (Fig. 1 and Ta- 
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blel) indicates from the peak areas that the b 
band should have ~1/10 intensity of the 5450 A 
band. The intensity in the gas phase is distributed 
between the three K components (Fig. 3). It has 
to be decided what is the appropriate rotational 
temperature. As is seen in Fig. 3, at 5 K the mid- 
dle peak, A"=l<— component, will be the most 
intense; at 20 K the outer two (K=2<—1 and 0<— 1) 
already become stronger. This will thus influ- 
ence a comparison of the DIB EW with laboratory 
data. 

The column density of Z-C3H2, ~5x 10 14 cm -2 — 
inferred in the previous section, was for the 5450 A 
DIB towards HD 183143. For the observation to- 
wards HD 206267 (Section 4.1) the EW=0.13A is 
lower and hence A(/-C3H 2 )«2x 10 14 cm~ 2 . We 
can now use this value to estimate the EW of 
the 6250 A band system (Fig. 3). The oscilla- 
tor strength of b is taken as /=0.0003, guessed 
by comparison with the corresponding absorption 
bands in the neon matrix. An uncertainty arises 
because of variation in phonon broadening and 
site structure (Fig. 1). Then at 6250 A this leads 
to an EW of around 10 mA. Thus below 5 K only 
the A=l<— peak would be seen and have EW of 
this order of magnitude. At a higher temperature, 
10— 30 K, this would be distributed between the 
three components K=2<- 1, A=l<— 0, K=0<- 1, 
each with EW of a few mA. 

An interesting inference concerns the anti- 
correlation with C3 abundance in diffuse clouds. 
From our measurements of the 4881 and 5450 A 
DIBs in four stars covering a large range in red- 
dening and supplemented by measurements of 



4881 and 5450 A DIBs. This is in contrast to the 
family of enhanced DIBs found by |Thorburn et| 
all d2003b for HD 204827 which correlates closely 



HD 204827 fromlHobbs et al. (2008), we find that 



these two DIBs correlate much more closely with 
each other than with E(B — V). This implies 
a common carrier. iHobbs et all (2008) specifi- 



cally observed HD 204827 because of its excep- 
tionally high column density of interstellar C3. 
They could detect none towards HD 183143. By 
contrast, from their spectra, the 4881 and 5450 A 
DIBs in HD 204827 are less than one third their 
strength in HD 183143. There is a similar, nega- 
tive trend between the line of sight abundance of 
interstellar C3 ( Oka et al|2003 ) and the strengths 
of the 4881 and 5450 A DIBs for the stars in our 
sample. The inference seems clear, hydrogenation 
of C3 to species such as /-C3H2 leads to deple- 
tion of C3 with a consequent enhancement of the 



with A(C 2 ). 

The present study provides laboratory and as- 
tronomical evidence for Z-C3H2 in the diffuse in- 
terstellar clouds and its identification as a DIB 
carrier. 
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4500 5000 5500 6000 

Wavelength (A) 



Fig. 1. — Absorption spectra of Z-C3H2 measured 
in 6 K neon matrices. The top trace was obtained 
by 1470 A photolysis of allene embedded in the 
matrix. It is identical with the published spec- 
trum in (|Hodges et al. 2000) generated by an- 



other chemical pathway. The middle and bottom 
traces were observed using mass-selected deposi- 
tion as described in the experimental section. La- 
bels a, b,c refer to the rotationally resolved gas- 
phase measurements ( A^feso^eFS]2006 Birza 
et~aT|[2005l >. 




5480 5450 5470 5490 4840 4880 4920 




5430 5450 5470 5490 A 4840 4880 4920 A 



Fig. 2. — Gas-phase laboratory spectra recorded 
by cavity ring-down spectroscopy in a supersonic 
slit jet expansion discharge using allene (traces 
a, c, d) and DC2D (trace b) as precursors. The 
continuous, solid lines are Gaussians fitted to the 
data, after partial removal of the C2 lines and sub- 
traction of the known Z-C3H absorptions (trace c). 
Absorptions in red are due to Z-C 3 H 2 . 
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6235 6350 6265 A 6150 6160 6170 A 

Fig. 3. — Simulated rotational profiles of the b 
(left) and c bands (right) of ^-C3ri2 which have 
been observed in the gas phase earlier. The struc- 
ture is characteristic of an asymmetric top, show- 
ing the three K components 2^—1, 1^—0 and 0^—1. 
The plot on the left uses 0.15 A linewidth for the 
individual rotational transition, corresponding to 
40,000 resolving power used in the astronomical 
observations. On the right 0.4 A is used, compat- 
ible with the conditions of the astronomical ob- 
servations as the pattern and broadening of the 
potassium atomic lines show. 
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5420 5440 5460 5480 A 



Fig. 4. — Simulation of the 5450 A band profile of 
Z-C 3 H 2 B 1 B 1 <-X 1 A 1 transition at 10 and 60 K 
taking 1 or 10 A linewidth. The laboratory absorp- 
tions shown in Fig. 2 were measured with a rota- 
tional temperature of around 60 K, determined by 
the rotational line intensity of the I-C3H2 bands 
above 6000 A. The 10 A linewidth corresponds to 
the 100 fs lifetime of the excited electronic state. 
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Fig. 5.— The 4881 A DIB, already known from the 
literature, in spectra of HD 46711, 190603, 183143 
and 7 0ri (unreddened standard), is compared to 
a Gaussian fit (FWHM 25 A) to the laboratory 
absorption spectrum of Z-C3H2 centred at 4887 A 
(dashed, red line). The short wavelength side of 
the DIB profile is distorted to varying degrees by 
H/3 at 4861 A. 
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5440 5450 5460 

Wavelength (A) 



Fig. 6.— Spectra of HD 46711, 206267, 186745, 
183143 and AE Aurigae showing the broad DIB at 
5450 A. The scales for the lightly reddened stars 
HD 206267 and AE Aurigae have been expanded 
by a factor 4. The gas-phase lab spectrum of l- 
C3H2 is shown by the dashed, red curve superim- 
posed on the HD 183143 spectrum. This confirms 
the agreement found earlier (Linnartz et al.|2010). 
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Fig. 7.— Interstellar KI line (7699 A) profiles 
for the reddened target stars aligned on their 
strongest interstellar component. The KI radial 
velocities actually adopted are given in Table 3. 
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Fig. 8.— Spectra of HD 204827, C Oph, the unred- 
dened standard, a Leo and model line profiles for 
Z-C 3 H 2 at 10 K and 0.15 A FWHM. The spectrum 
for £ Oph has been expanded by a factor of 5 rel- 
ative to the others. The broad feature centred at 
6159 A in a Leo is a rotationally broadened stellar 
line. Narrow DIB lie within the 6159 and 6166 A 
model profiles for the two reddened stars. 
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Fig. 9.— Spectra of HD 204827, HD 206267, 
COph, and HD 179406 together with model line 
profiles for Z-C 3 H 2 at 10 K and 0.15 A FWHM. The 
spectrum for £ Oph has been expanded by a factor 
of 2.5 while that for HD 204827 has been halved. 
There are DIB consistently within the limits of the 
model line profiles. 
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Table 1 

Absorption band maxima (A) of ;-C 3 H 2 and Z-C 3 D 2 measured in a 6K neon matrix and in the 

GAS PHASE (if=l«-0). 



Neon Matrix Gas Phase Transition Band 



A(J-C 3 H 2 ) 


1 


A((-C 3 D 2 ) 


I 


a(;-c 3 h 2 ) 


A(Z-C 3 D 2 ) 








6284 


0.2 


6260 


0.4 


6318.9 




A X A 2 


•i-A^Ai 


a 


6219 


0.3 


6130 


1.4 


6251.7 






II 


b 


6122 


0.9 


6093 


1.5 


6159.2 






II 


c 


5445/5417 1 


10 


5412 


10 


5450(3) 


5458(3) 


B x Bi 


4-X x Ai 


2 1 


5143 


3 


5160 


3.5 


5165-5185 2 






n 


2 1 4 1 


4856 


9 


4857 


12 


4887(3) 






n 


2 2 

„o 


4633 


4.4 






4645-4665 2 






ii 


2 2 4 1 


4412 


4.7 






4425-4445 2 






ii 


2 3 



most intense site 

extrapolated values based on neon/gas shifts of the 2j band 
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Table 2 

4881 and 5450 A DIB central depths (%) and EW (mA). 



Star 


s P 


E(B-V) 


4881 


5450 








% 


EW 1 


% 


EW 2 


HD 46711 


B3 II 


1.05 


5.7 


1425 


3.0 


390 


183143 


B7 la 


1.28 


6.6 


1650 


3.0 


390 


186745 


B8 la 


0.78 






2.6 


338 


190603 


B1.5 Ia+ 


0.70 


4.3 


1075 






AEAur 


09 V 


0.48 






0.8 


104 


206267 


06 


0.53 






1.0 


130 



1 based on the measured % central depth and a FWHM of 25 A 

2 based on the measured % central depth and a FWHM of 13 A 
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Table 3 
6000 a triplet targets. 



Star 


s P 


E(B - V) 


IS KI 1 


Spectrograph 








km s~ 1 


6150 


6250 


COph 


09.5 V 


0.30 


-14.53 


Esp 2 


Gee 3 


HD 179406 


B3 V 


0.27 


-12.53 




Esp 


204827 


BO V 


1.11 


-17.0 


Esp 


Esp 


206267 


06 e 


0.53 


-17.0 




Gee 



1 for strongest component — see Fig. 8 
2 ESPaDONs i?=80,000 
3 Gecko R=120,000 
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Table 4 

Comparison of laboratory and interstellar 6250 A triplet component wavelengths A (A) 1 
and EW (mA). The lab values give the maximum of the K=l<r- band using a Lorentzian 

WIDTH OF ~0.15 A. 



Source 


A 


A 


EW 


A 


EW 


A 


EW 


A 


EW 


A 


EW 


lab 


6152.8 


6159.0 




6165.8 




6244.4 




6251.6 




6259.0 




HD 179406 












6244.4 
6245.3 


2.4 
1.5 


6251.0 
6252.5 


1.6 
1.0 


6259.5 


2.2 


HD 206267 












6244.5 
6245.3 


2.3 
1.5 


6251.0 
6252.2 


1.2 
0.4 


6259.5 


1.1 


HD 204827 




6158.5 


8.0 


6165.8 


3.0 


6244.6 
6245.2 


5.0 
1.5 


6251.0 
6252.2 


4.3 
2.6 


6259.4 


3.0 


COph 




6158.5 


1.4 


6165.8 


1.6 


6244.5 
6245.3 


0.5 
0.3 


6251.0 
6252.2 


0.6 
0.3 


6258.4 
6259.6 


0.2 
0.5 



'some interstellar components are paired 
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